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HIGHLIGHTS 


►  A  new  method  is  proposed  to  describe  the  current— voltage  and  power— voltage  characteristics  of  solar  cell  and  module  in  this  article. 

►  By  using  the  exact  explicit  analytical  solutions,  the  I-V  curves  of  the  corresponding  modules  are  calculated  and  simulated. 

►  To  extract  the  parameters,  an  optimized  technique  is  presented  using  Lambert  W  function  and  polynomial  curve  fitting. 
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A  new  method  is  proposed  to  describe  the  characteristics  of  solar  cell  and  module  in  this  article.  By  using 
the  exact  explicit  analytical  solutions,  the  current-voltage  and  power-voltage  curves  of  the  corre¬ 
sponding  modules  are  calculated  and  simulated.  In  order  to  extract  the  parameters  of  the  solar  cells,  an 
optimized  technique  is  presented  by  using  Lambert  W  function  and  polynomial  curve  fitting.  The 
accuracy  of  the  model  is  compared  with  previous  methods  in  other  work.  Results  have  the  good 
agreement  between  the  fitted  current— voltage  curves  and  the  experimental  data.  Moreover,  error  and 
statistical  analyses  are  carried  out  to  illustrate  the  accuracy  of  the  estimated  parameters  and  procedure 
suitability. 

©  2012  Published  by  Elsevier  B.V. 


1.  Introduction 

Due  to  the  environmental  and  economic  benefits  obtained  from 
photovoltaic  generation,  many  research  works  have  addressed  the 
development  of  solar  power  system  in  recent  years.  One  of  the 
fundamental  issues  is  to  have  access  to  model  and  parameters  extrac¬ 
tion  for  the  solar  cells.  Over  the  years,  several  models  were  developed  to 
describe  the  non-liner  characteristics  of  the  current— voltage  (/-V) 
curve  of  the  solar  cells.  The  most  popular  approach  is  to  utilize  the 
electrical  equivalent  circuit,  which  is  primarily  based  on  the  Shockley 
diode  equation.  In  practical  applications,  the  single  diode  model  is  the 
main  equivalent  circuit  [1—6].  Moreover,  the  single  diode  model 
contains  five  parameters,  which  are  called  photocurrent,  saturation 
current,  series  resistance,  shunt  resistance  and  ideality  factor. 

However,  due  to  the  transcendental  nature  of  the  current  equation 
for  solar  model,  significant  computation  effort  is  required  to  obtain  all 
the  model  parameters  [7].  To  address  this  issue,  several  parameters 
extraction  methods  have  been  developed.  One  of  the  methods  is 
numerical  method,  which  typically  uses  polynomials  to  demonstrate 
the  current- voltage  relationship  of  PV  panels  and  implements  iterative 


*  Corresponding  author.  Tel./fax:  +86  021  67792587. 

E-mail  address:  sunyz@dhu.edu.cn  (Y.  Sun). 

0378-7753 /$  -  see  front  matter  ©  2012  Published  by  Elsevier  B.V. 
http://dx.doi.Org/10.1016/j.jpowsour.2012.07.061 


method  such  as  the  recursive  least-squares  method  [7]  and 
Newton-Raphson  method  (NRM)  [8]  to  obtain  all  the  model  param¬ 
eters.  However,  the  disadvantage  of  this  approach  is  its  dependency  on 
the  initial  values  used  in  the  proposed  iterative  technique.  Moreover, 
NRM  does  not  always  converge.  Another  solution  technique  is  called 
analytic  method  [9]  which  is  presented  to  express  the  transcendental 
current-voltage  characteristic  containing  parasitic  power  consuming 
parameters  like  series  and  shunt  resistances.  Lambert  W  function  is 
always  applied  to  solve  the  transcendental  equations  [7,9-13],  and  has 
been  widely  used  in  many  fields.  For  example,  R.M.  Corless  et  al.  have 
given  many  available  applications  on  the  Lambert  W  function  [14].  H. 
Abebe  et  al.  used  Lambert  W  function  to  develop  an  analytical  compact 
model  for  the  asymmetric  lightly  doped  Double  Gate  (DG)  MOSFET  [9]. 
However,  fewer  researchers  have  given  all  the  analytical  solution  for 
the  five  parameters  of  the  single  diode  model.  Although  the  analytical 
resolution  is  given,  the  data  is  still  less  accurate  and  common.  Recently, 
various  high  accuracy  algorithms  techniques  have  been  reported,  such 
as  particle  swarm  optimization  (PSO)  [15],  differential  evolution  (DE) 
[16],  genetic  algorithm  (GA)  [17]  and  pattern  search  (PS)  [18,19].  While 
the  calculation  process  is  more  complicated,  this  is  not  conducive  to 
master. 

Considering  the  many  advantages  of  numerical  and  analytic 
method,  this  paper  proposes  a  new  method  for  determining  the  I-V 
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Nomenclature 

T 

cell  temperature  (I<) 

PMAE 

maximum  absolute  power  error  (W) 

CMAE 

maximum  absolute  current  error  (A) 

PMAEmp 

maximum  absolute  power  error  at  the  maximum 

CMAEmp 

maximum  absolute  current  error  at  the  maximum 

power  point  (W) 

power  point  (A) 

PRMSE 

root  mean  squared  error  at  the  maximum  power  point 

CRMSE 

root  mean  squared  error  at  the  maximum  power  point 

(W) 

(A) 

Q 

electron  charge  (1.6  x  10  19C) 

cP 

the  conductance  of  the  shunt  resistance(Q  1) 

V 

output  voltage  of  the  cell  (V) 

/ 

output  current  of  the  cell  (A) 

VCi 

computed  values  of  the  I—V  curves  for  a  given  voltage 

II 

photocurrent  (A) 

point 

fmp 

current  at  the  maximum  power  point  (A) 

Vei 

experimental  values  of  the  I—V  curves  for  a  given 

he 

short  circuit  current  (A) 

voltage  point 

Io 

saturation  current  of  the  equivalent  diode  (A) 

Vmp 

voltage  at  the  maximum  power  point  (V) 

k 

Boltzmann  constant  (1.38  x  10  23  J  K-1) 

Vac 

open  circuit  voltage  (V) 

n 

ideality  factor  of  the  cell 

6 

angle  between  Voc  and  Jsc  (o) 

N 

number  of  measurements 

0 1 

angle  between  Voc  and  Vmp  (o) 

Rs 

series  resistance  of  the  cell  (Q) 

@2 

angle  between  /mp  and  Jsc  (o) 

Rp 

shunt  resistance  of  the  cell  (Q) 

characteristics  of  solar  cells.  To  validate  the  accuracy  of  the  new 
method,  sets  of  comparative  results  were  tested  in  matlab/simulink 
environment.  Comparing  with  other  competing  methods,  the  ob¬ 
tained  results  have  high  accuracy.  Meanwhile,  the  comparative 
results  show  that  the  proposed  method  gives  a  good  representation 
at  the  maximum  point. 


2.  Theory  and  application 

2.1.  Explicit  analytic  solutions 

For  single  diode  model  of  solar  cells,  the  I-V  characteristics  can 
be  expressed  as: 


F(I,  V)  =  j  V  d/ 

0 

=  T[(_V-Bj  +  c)2-(-V  +  c)2]  -1 BI2+ADI  (6) 

Here  the  parameters  of  A ,  B,  C,  D  are  given  by  the  following 
expressions: 


RP 

(7) 

Rs  +  Fp 

(8) 

nl/th  +  Fp(4  + 10) 

(9) 

I  =  Il-Io 


exp 


V  +  IR< 
nVlh 


1 


V  +  IRs 

Ro 


(1) 


Here  Vih  is  the  “thermal  voltage”  of  a  cell,  which  is  given  by: 
Vth(T)  =  kT/q  (2) 


D  =  /L  +  /o 

From  Eq.  (1)  and  Eqs.  (7)-(10),  it  is  obtained: 
RP  =  A 


(10) 


(11) 


The  solution  of  Eq.  (1)  can  be  expressed  by  using  Lambert  W  Rs  =  B  -  A 
function  [10]: 


(12) 


/  = 


Rpik+M-V 


Rs  +  Rp 

.HlW _ 

Rs  \nVth(R, 


RsRplo 


Rr 


-exp 


RpjRA  +  Rslo  +  V) 


nKh(Rs  +  Rp) 


(3) 


v  -  '>  -  *•'  -  } 

(4) 

Similarly,  the  output  power  can  be  given  explicitly  in  terms  of  V 
as  follows: 


P(V) 


m/l+w-v2 

Rs  “F  Rp 

_  nWth  f  RsRpIo  Rp(Rsk  +  Rslo  +  V) 

Rs  lnVth(^s  +  FP)  _  ft^th(Fs  +  Fp) 


(5) 


The  integration  of  Eq.  (4)  can  be  expressed  as  below  [20]: 


n  =  ( C-AD)/Vth 


D  —  I  — 


V  +  I(B-A ) 


In  = 


exp 


V  +  I(B-A) 


C-AD 


-2 


(13) 


(14) 


4  =  D-Io  (15) 

Once  A,  B,  C,  D  are  determined,  the  five  parameters  that  describe 
solar  cells  model  behavior  are  calculated  in  sequence  with  Eqs. 
(11)— (15),  respectively.  Once,  the  five  parameters  are  given, 
through  Eqs.  (3)  and  (5),  the  I-V  and  P—V  curves  of  the  corre¬ 
sponding  modules  are  calculated  and  simulated. 


2.2.  Polynomial  curve  fitting 

Eq.  (6)  shows  that  the  integral  of  the  solar  voltage  F(J,V)  is  the 
polynomial  function  on  A ,  B,  C  and  D.  And  the  F(/,V)  can  be  obtained 
from  the  experimental  data  by  polynomial  curve  fitting.  Therefore, 
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substituting  four  different  groups  of  value  F(J,V)  into  Eq.  (6),  the  A,  B, 
C  and  D  can  be  determined  by  solving  the  algebraic  equations. 

In  order  to  accurately  represent  the  characteristics  of  solar  cells, 
a  sixth-order  polynomial  [8]  is  defined  as: 

V(Ii)  =  a0  +  fll  +  a2lf  +  a3lf  +  a4l?  +  a5lf  +  a6lf  (16) 


seven  experimental  data  points  which  were  used  for  polynomial 
curve  fitting,  can  be  chosen  by  the  following  formula,  as  shown  in 
Fig.  1. 

J0  =  tan^/mp 


Substituting  Eq.  (16)  in  Eq.  (6),  F(J,V)  can  be  obtained: 

F(I,V)  =  a0l+l-ail2 +1  a2I3  +1  a3I4  +±a4I5 +±a5I6 +la6/7 

=  2^[(-V-B/+C)2-(-V+C)2]  ~BI2+ADI  (17) 


h  = 


mp 


1  -^(1  -  cos  6) 


h  —  Imp  ^  1  1  COS  ft )  Imp 


For  evaluating  the  error  of  this  method,  the  measure  of  error  is 
given  by: 


S  =  \\e\\i  = 


-  va2  =  X>tt)  -  va2 


\  1  =  0 


(18) 


1  =  0 


To  minimize  the  error,  the  derivative  with  respect  to  the  sample 
points  is  zero. 

=  =  0,  ./  =  0.1 . 6  (19) 


i= 0 \k= 0 


And  Eq.  (19)  can  be  written  as: 

6/6  \  m 

E  E4+fck  =  E^>  j  =  e 


i=0 \k= 0 


(20) 


1  =  0 


For  convenience  of  representation,  the  following  notations  of 
vectors  and  matrix  are  introduced: 


Y  = 


6 

E  vi 

i  =  0 

E  Wi 

i  =  0 
6  ‘ 

E  if  Vf 

Li  =  0 


X  = 


6  +  1 

Ei,  - 

i-0 

•  Elf 

i  =  0 

6 

E  h 

1=0 

6 

E'f  ■■ 
1  =  0 

6 

•  EiJ 

1  =  0 

6'  _ 
E  If 

Li  =  o 

Ef  ■ 

i  =  0 

•  Ei?2 
1  =  0 

~a0~ 

ai 

_fl6_ 

6  = 


Then,  the  measured  output  vector  y  can  be  represented  by 
a  simple  regression  model: 


y  =  xe 


(21) 


According  to  the  theorem  of  the  matrix  transformation,  the 
vector  of  the  estimated  parameters  6  is  given  by: 


6  =  X-’Y 


(22) 


Due  to  the  existing  shape  similarity  between  the  characteristics 
of  solar  cells  and  the  step  response  of  a  first  order  system  [23],  the 


h  =  h 


mp 


J4  —  Imp  +  "7(1—  COS  0  )  (Jsc  —  /mp) 


,  _  1,  1 

h  ~  2hc~2 


fmp  +  (1  —  COS  6)  (Jsc  ~  I] 


mPj 


I6  =  tan  d2hi 


Previous  studies  have  shown  that  Rs  impacts  the  I—V  curve  near 
the  maximum  power  point,  and  Rp  determines  the  slope  of  I—V 
curve  between  V0c  and  Vmp.  A  typical  approach  [5]  is  to  estimate 
Rs  and  Rp  value  by  using  the  slopes  at  the  Voc  and  /sc,  respectively.  Jo 
and  JL  are  always  calculated  at  the  points  of  (V0o  0)and  (0,JSC).  And  n 
is  determined  by  the  inherent  characteristics  of  the  solar  cells  [2]. 
Thus,  to  accurately  obtain  the  five  parameters,  the  four  different 
groups  of  value  F(J,V)  are  chosen  at  Jo,  h ,  UJe-  Substituting  the  four 
groups  of  value  F(J,V)  into  Eq.  (17),  the  A,  B ,  C  and  D  can  be  deter¬ 
mined  by  solving  the  algebraic  equations. 


3.  Validation  and  discussions 

In  an  attempt  to  confirm  the  validity  of  our  method,  a  practical 
data  of  a  silicon  solar  cell  and  module  was  used.  And  the  experi¬ 
mental  data  was  reported  by  Easwarakhanthan  et  al.  [21].  In 
addition,  another  different  type  of  solar  cell  called  plastic  solar  cell 
was  used  [22].  Moreover,  the  proposed  algorithm  was  implanted  in 
Matlab/Simulink  environment.  Sixth-order  polynomial  curve 
fitting  and  well-picked  four  groups  of  value  F(J,V)  were  employed 
for  parameters  extraction  from  commercial  solar  cells.  To  provide 
a  comprehensive  evaluation,  the  obtained  results  were  also 
compared  with  other  methods  by  using  the  same  reference 
parameters,  and  they  were  given  in  Tables  1-3.  Through  Eqs.  (3) 


Fig.  1.  The  seven  experimental  data  used  for  polynomial  curve  fitting. 
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Table  1 

Parameters  for  the  silicon  solar  cell  and  module  using  different  methods. 


Item 

Silicon  solar  cell  (33  °C) 

Silicon  solar  module  (45  °C) 

Ref.  [19] 

Ref.  [20] 

Ref.  [16] 

Our  method 

Ref.  [19] 

Ref.  [20] 

Ref.  [16] 

Our  method 

4(A) 

0.7608 

0.7603 

0.7617 

0.7609 

1.0318 

1.0300 

1.0313 

1.0313 

h  (fiA) 

0.3223 

0.3374 

0.9980 

0.3220 

3.2876 

6.3986 

3.1756 

3.2212 

n 

1.4837 

1.4841 

1.6000 

1.4837 

48.4500 

50.9900 

48.2889 

48.3221 

Rs  (A) 

0.0364 

0.0376 

0.0313 

0.0364 

1.2057 

1.1619 

1.2053 

1.2132 

Cp(0-’) 

0.0186 

0.0094 

0.0156 

0.0185 

0.0018 

0.0014 

0.0014 

0.0016 

Table  2 

Statistical  results  for  the  silicon  solar  cell  and  module  using  different  methods. 

Item 

Silicon  solar  cell  (33  °C) 

Silicon  solar  module  (45  °C) 

Ref.  [20] 

Ref.  [16] 

Our  method 

Ref.  [20] 

Ref.  [16] 

Our  method 

CMAE 

0.01550 

0.02819 

0.01258 

0.01731 

0.004423 

0.005152 

CMAEmp 

6.998  x  10"5 

7.934  x  10"3 

1.551  x  10"3 

0.0090 

0.0043 

0.0007 

CRMSE 

0.0051 

0.0091 

3.543  x  10"3 

0.0314 

0.0028 

0.0026 

PMAE 

0.008834 

0.01607 

0.007174 

0.2770 

0.05394 

0.08244 

PMAEmp 

3.149  x  10"5 

3.570  x  10"3 

6.980  x  10"3 

0.1098 

0.0535 

0.0089 

PRMSE 

0.0026 

0.0048 

9.552  x  10"4 

0.3865 

0.0325 

0.0382 

Table  3 

Parameters  and  statistical  results  for  the  plastic  solar  cell  using  different  methods. 

Plastic  solar  cell(27.3  °C) 

Item 

/l(A)  /o(mA) 

n 

Rs(fi)  Gp(cr’) 

CMAE  (pA)  CMAEmp 

(pA)  CRMSE  (pA) 

PMAE  ( pW)  PMAEmp  ( pW) 

PRMSE  (pW) 

Ref.  [13] 

0.0079  0.0136 

2.31 

8.59  0.0051 

0  0 

0 

0 

0 

0 

Ref.  [22] 

0.0079  0.0329 

2.59 

8.586  0.0050 

920.93  106.74 

288.89 

841 

59.775 

205.89 

Our  method 

0.0079  0.0136 

2.3101 

8.5884  0.0051 

1.33  0.188 

0.394 

0.999 

0.106 

0.277 

and  (5),  the  I—V  and  P—V  curves  of  the  corresponding  modules  are 
calculated  and  simulated.  The  results  for  both  the  silicon  solar  cell 
and  module  are  given  in  Figs.  2  and  3.  The  errors  in  different 
methods  are  shown  in  Figs.  4  and  5.  Moreover,  the  results  for  the 
plastic  solar  cell  are  given  in  Fig.  6. 

With  Tables  1  and  3,  it  is  easy  to  see  that  parameters  extracted  by 
using  the  proposed  method,  are  very  close  to  those  reported  in  the 
other  two  references.  Figs.  2, 3  and  6  show  that  the  new  method  can 
achieve  high  accuracy,  which  are  in  good  agreement  with  the 
experimental  results.  In  order  to  make  this  study  as  general  as 
possible,  the  same  statistical  analysis  is  performed  [17].  The  absolute 
error  (AEj)  and  the  root  mean  squared  error  (RMSE)  are  given  by: 

AEj  =  \Vti  -  Veil  (17a) 


RMSE  = 


\ 


i  N 


1=1 


(18a) 


Figs.  3  and  4  show  that,  the  new  method  can  represent  the  cell 
and  module  characteristics  with  the  small  absolute  errors  of 
current  and  power.  In  addition,  it  is  interesting  to  find  that,  in 
different  voltage  ranges  the  absolute  error  varies  with  different 
methods.  Almost  before  the  maximum  power  point,  the  error  of 
three  algorithms  is  basically  the  same,  after  which  there  is  a  clear 
distinction.  Moreover,  comparisons  with  simulation  results,  it  has 
shown  that  the  proposed  method  has  a  smaller  error  value  after  the 
maximum  power  point.  The  results  indicate  that  the  parallel 
resistance  has  a  greater  impact  on  the  maximum  power  point. 

In  addition,  it  can  be  concluded  from  Table  2  that,  this  new 
method  can  accurately  access  the  I-V  characteristics  of  the  solar 
cell  and  module  with  much  lower  errors.  For  the  maximum  voltage 
and  current,  CMAEmp  and  PMAEmp  values  are  less  than  0.0043  A 
and  0.1098  W,  respectively.  The  root  means  squared  errors  at 
current  and  power  do  not  exceed  0.01258  A  and  0.0382  W. 
Comparisons  with  simulation  results  under  different  methods  in 
Table  2,  it  has  shown  that  the  proposed  method  reduces  the  overall 
estimation  error  at  a  very  low  standard  deviation. 


Fig.  2.  Experimental  data  and  the  fitted  curve  for  silicon  solar  cell:  (a)  I-V  characters  and  (b)  P—V  characters. 
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Experimental  data  — 
Fitted  curve 


Fig.  3.  Experimental  data  and  the  fitted  curve  for  silicon  solar  module:  (a)  I-V  characters  and  (b)  P—V  characters. 


Voltage  (V) 

Fig.  4.  Absolute  errors  for  silicon  solar  cell:  (a)  absolute  current  errors  and  (b)  absolute  power  errors. 


Curve  for  reference[20]- 
Curvc  for  rcfcrcncc[  1 6]  - 
Fitted  curve  - 


|  Curve  for  refcrencc[20] 

w  00i(-  Curve  for  rcfcrcnccl  16] 


Fig.  5.  Absolute  errors  for  silicon  solar  module:  (a)  absolute  current  errors  and  (b)  absolute  power  errors. 


Experimental  data  - 
Fitted  curve  ■ 


01  02  03  04  0$ 

Voltage  (V) 


00  07  OB 
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Experimental  data  -  _»«»** 
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■ 

\ 

\ 

\ 

.  .  \  1 

Fig.  6.  Experimental  data  and  the  fitted  curve  for  plastic  solar  cell:  (a)  I—V  characters  and  (b)  P—V  characters. 
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Therefore,  there  are  three  aspects  which  have  been  addressed. 
The  first  question  involves  the  I—V  characteristics,  the  second 
problem  relates  to  the  parameters  extraction,  the  third  aspect  deals 
with  the  error  and  statistical  analyses.  We  tentatively  put  forward 
that,  the  quality  standard  of  a  good  method  has  not  only  a  small 
current  error,  but  also  small  deviation  of  power,  especially  at  the 
maximum  power  point. 

4.  Conclusions 

We  have  presented  a  new  method  for  determining  the  I—V  and 
P—V  characteristics  of  solar  cells.  By  using  the  exact  explicit  analyt¬ 
ical  solutions,  the  I—V  and  P—V  curves  of  the  corresponding  modules 
are  calculated  and  simulated.  A  new  technique  is  presented  by  using 
Lambert  W  function  and  subsection  polynomial  curve  fitting  to 
extract  the  parameters.  The  method  has  been  successfully  applied  to 
commercial  solar  cell  and  module.  The  results  obtained  are  in  good 
agreement  with  those  published  previously.  In  addition,  the  quality 
standard  of  a  good  method  has  not  only  a  small  current  error,  but  also 
small  deviation  of  power,  especially  at  the  maximum  power  point. 
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